Filling dependence of the Mott transition in the degenerate Hubbard model 
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Describing the doped Fullerenes using a generalized Hubbard model, we study the Mott transition 
for different integer fillings of the ti u band. We use the opening of the energy-gap E g as a criterion 
for the transition. E g is calculated as a function of the on-site Coulomb interaction U using fixed- 
node diffusion Monte Carlo. We find that for systems with doping away from half-filling the Mott 
transitions occurs at smaller U than for the half-filled system. We give a simple model for the doping 
dependence of the Mott transition. 
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I. INTRODUCTION 



The Hubbard Hamiltonian is a simple model for study- 
ing strongly interacting systems. In particular it is used 
to investigate the Mott-Hubbard metal-insulator transi- 
tion in half-filled systems^ It is clear that for strong cor- 
relations such a system should be insulating, since in the 
atomic limit the states with exactly one electron per lat- 
tice site are energetically favored, while all other states 
are separated from those by a Coulomb gap. For a gener- 
alized Hubbard model with degenerate orbitals the same 
argument implies that for strong correlations not only 
the half-filled, but all integer filled systems will become 
Mott-Hubbard insulators. It is then natural to ask how 
the location of the transition depends on the filling. 

As an example we consider a Hamiltonian describing 
the alkali doped Fullerides. □ It comprises the three-fold 
degenerate t\ u orbital and the Coulomb interaction U 
between the electrons on the same molecule. Using this 
Hamiltonian, we have recently shown that, although U 
is substantially larger than the band width W, K3C60 
not a Mott insulator but a (strongly correlated) metal.El 
Prompted by the synthesis of an isostructural family of 
doped Fullerenes A„C6o with different fillings n,u we now 
address the question of the Mott transition in integer 
doped Fullerides. For these systems we have the inter- 
esting situation that for fillings n = 1, 2, 3, 4, and 5 cal- 
culations in the lo£al density approximation predict them 
all to be metallic Jj while in Hartree-Fock they all are in- 
sulators. Performing quantum Monte Carlo calculations 
for the degenerate Hubbard model at different fillings and 
for values of U typical for the Fullerides, we find that all 
the systems are close to a Mott transition, with the crit- 
ical correlation strength U c at which the transition takes 
place strongly depending on the filling n. More generally, 
our results show how, for an otherwise identical Hamil- 
tonian, the location of the Mott transition U c depends 
on the filling. U c is largest at half-filling and decreases 
for fillings smaller or larger than half. We contrast these 
findings with the results from Hartree-Fock calculations 
which predict a much too small U c and show almost no 



doping dependence. We give an interpretation of the re- 
sults of the quantum Monte Carlo calculations extending 
the hopping argument introduced in Ref. [| to arbitrary 
integer fillings. Despite the crudeness of the argument it 
explains the doping dependence found in quantum Monte 
Carlo. We therefore believe that our simple hopping ar- 
gument captures the basic physics of the doping depen- 
dence of the Mott transition in degenerate systems. 

In section II we introduce the model Hamiltonian for 
doped Fullerenes with a three-fold degenerate t\ u band. 
We discuss the fixed-node approximation used in the dif- 
fusion Monte Carlo calculations, present the results of our 
quantum Monte Carlo calculations, and contrast them 
to the result of Hartree-Fock calculations. Section III 
gives an interpretation of the results of our calculations 
in terms of intuitive hopping arguments. We introduce 
the many-body enhancement of the hopping matrix el- 
ements, which explains how orbital degeneracy N helps 
to increase the critical U at which the Mott transition 
takes place and we analyze how frustration leads to an 
asymmetry of the critical U for fillings n and 2N — n. A 
summary in Sec. IV closes the presentation. 



II. MODEL CALCULATIONS 



A. Model Hamiltonian 



Solid Cgo is characterized by a very weak inter- 
molecular interaction. Therefore the molecular -levels 
merely broaden into narrow, well separated bands B The 
conduction band originates from the lowest unoccupied 
molecular orbital, the 3-fold degenerate ti u orbital. To 
get a realistic, yet simple description of the electrons in 
the ti u band, we use a Hubbard-like model that describes 
the interplay between the hoppuag of the electrons and 
their mutual Coulomb repulsion!] 
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The sum (ij) is over nearest-neighbor sites of an fee lat- 
tice. The hopping matrix elements ti m j m i between or- 
bital m on molecule i and orbital ml on molecule-,j are 
obtained from a tight-binding parameterization!] The 
molecules are orientationally disordered,!!] and the hop- 
ping integrals are .chosen such that this orientational dis- 
order is included. □ The band width for the infinite sys- 
tem is W — 0.63 eV. The on-site Coulomb interaction is 
U 1=3 1.2 eV . The model neglects multiplet effects, but we 
remark that these tend to be counteracted by the Jahn- 
Teller effect, which is also not included in the model. 

We will investigate the above Hamiltonian for differ- 
ent integer fillings n of the t\ u band. The correspond- 
ing Hamiltonians describe a hypothetical family of doped 
Fullerides A„Ceo with space group Fm3m, i.e. an fee lat- 
tice with orientationally disordered Ceo molecules. In 
the calculations we use the on-site Coulomb interaction 
U as a parameter to drive the system across the Mott 
transition. 



B. Quantum Monte Carlo method 

As the criterion for determining the metal-insulator 
transition we use the opening of the gap 

E g = E(N + l)-2E(N) + E{N-l), (2) 

where E(N) denotes the total energy of a cluster of N mo i 
molecules with N electrons in the t\ u band. Since we are 
interested in integer filled systems, N = nN mo i, n an 
integer. For calculating the energy gap (|^) we then have 
to determine ground-state energies for the Hamiltonian 
(|l|). This is done using quantum Monte Carlo.EI Starting 
from a trial function \^>t) we calculate 

\yW} = [l-T(H- W )} n \y T ), (3) 

where w is an estimate of the ground-state energy. The 
|^(")) are guaranteed to converge to the ground state 
l^o) of H, if r is sufficiently small and \^t) is not or- 
thogonal to l^o)- Since we are dealing with Fermions, 
the Monte Carlo realization of the projection (|^) suf- 
fers from the sign-problem. To avoid the exponential 
decay of the signal-to-noise ratio we use the fixed-node 
approximations For lattice models this involves defining 
an effective Hamiltonian H e g by deleting from H all non- 
diagonal terms that would introduce a sign-flip. Thus, by 
construction, H g is free of the sign-problem. To ensure 
that the ground-state energy of H c g is an upper bound 
of the ground state of the original Hamiltonian H, for 
each deleted hopping term, an on-site energy is added in 
the diagonal of H c g. Since Y$>t) is used for importance 
sampling, iJ c ff depends on the trial function. Thus, in a 
fixed-node diffusion Monte Carlo calculation for a lattice 



Hamiltonian, we choose a trial function and construct 
the corresponding effective Hamiltonian, for which the 
ground-state energy -Efndmc can then be determined by 
diffusion Monte Carlo without a sign-problem. 

For the trial function we make the Gutzwillcr Ansatz 

\nu ,g)) =g D muo)), (4) 

where the Gutzwiller factor reflects the Coulomb term 
U D = U ni ma m m io' in the Hamiltonian (p. |$(£/o)) 
is a Slater determinant that is constructed by solving 
the Hamiltonian in the Hartree-Fock approximation, re- 
placing U by a variational parameter Uq. Details on the 
character of such trial functions and the optimization of 
Gutzwiller parameters can be found in Ref. [l0| 

TABLE I. Total energy (in eV) for a cluster of four Ceo 
molecules with 6 + 6 electrons (filling 3) for different values of 
the on-site Coulomb interaction U. The difference between 
the fixed-node diffusion Monte Carlo results and the exact 
ground-state energy is shown in the last column. Note that 
Efndmc is always above the exact energy, as expected for a 
variational method. 



-Efndmc A.E 

025 0.8457 0.8458(1) 0.000 

0.50 4.1999 4.2004(1) 0.001 

0.75 7.4746 7.4756(2) 0.001 

1.00 10.6994 10.7004(2) 0.001 

1.25 13.8860 13.8875(3) 0.002 

1.50 17.0408 17.0427(4) 0.002 

1.75 20.1684 20.1711(5) 0.003 

2.00 23.2732 23.2757(10) 0.003 



TABLE II. Total energy (in eV) for a cluster of four Ceo 
molecules with on-site Coulomb interaction U — 1 eV for dif- 
ferent number of electrons iVf + iVj. . The difference between 
the fixed-node diffusion Monte Carlo results and the exact 
ground-state energy is shown in the last column. -Efndmc is 
always above the exact energy, as expected for a variational 
method. 
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8.4649 


8.4677(2) 


0.003 
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10.6994 


10.7004(2) 


0.001 
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13.3973 


13.3973(2) 


0.001 
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19.5094 


19.5109(3) 


0.002 
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22.9515 


22.9530(3) 


0.002 
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26.6590 


26.6613(3) 


0.002 



To check the accuracy of the fixed-node approximation, 
we have determined the exact ground-state energies for 
a (small) cluster of four Gqo molecules using the Lanczos 
method. For systems with different on-site Coulomb in- 
teraction (Table ||) and varying number of electrons (Ta- 
ble R) , we consistently find that the results of fixed-node 
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diffusion Monte Carlo are only a few meV above the ex- 
act energies. 



C. Quantum Monte Carlo results 



Since the quantum Monte Carlo calculations are for fi- 
nite clusters of N mo i molecules, we have to extrapolate 
the calculated energy gaps to infinite system size. An 
obvious finite-size effect is the fact that the one-particle 
spectrum is discrete, hence there can be a gap, even for 
U = 0. Furthermore, in evaluating (^), we add and sub- 
tract one electron to a finite system. Even if we distribute 
the extra charge uniformly over all molecules, there will 
be an electrostatic contribution of U/N mo i to the gap. 
We therefore introduce 



E G = E g - E g {U = 0) 



U 



N, 



mol 



(•5) 



These corrections are expected to improve the finite-size 
extrapolation. In practice they turn out to be quite small. 
For a cluster of 32 Cgo molecules, e.g., E g (U — 0) is typ- 
ically already less than 10 meV. In the thermodynamic 
limit both correction terms vanish, as they should. 

The results of the quantum Monte Carlo calculations 
are shown in Fig. [|. Plotting the finite-size corrected 
gap Eg for different values of the Coulomb interaction U 
versus the inverse system size l/N mo i, we read off where 
the gap starts to open. For the system with one electron 
per molecule the gap opens around U c « 0.75 . . . 1.00 eV. 
At filling 2 the transition takes place later, at U c « 
1.25... 1.50 eV. For both, filling 3 and 4 we find the 
largest critical U: U c « 1.50 ... 1.75 eV. For the sys- 
tem with 5 electrons per molecule the gap opens around 
U c « 1.00 . . . 1.25 eV. The results are summarized in Fig. 
||. Thus we find that for an otherwise identical Hamil- 
tonian the critical U for the Mott transition depends 
strongly on the filling. U c is largest at half-filling and 
decreases away from half-filling. The decrease in U c is, 
however, not symmetric around half-filling. It is more 
pronounced for fillings < 3 than for fillings > 3. 

We note that the opening of the gap is accompa- 
nied by a change in the character of that trial function 
which yields the lowest energy in the fixed-node approx- 
imation. For small [/, where the system is still in the 
metallic regime, paramagnetic trial functions with small 
Uq (see the discussion after eqn. (|])) are best. When 
the gap starts to open, trial functions with larger Uq, 
which have antiferromagnetic character, give lower ener- 
gies. The corresponding Slater determinants |<&(f/o)) de- 
scribe a Mott insulator in Hartree-Fock approximation. 
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FIG. 1. Finite-size corrected gap Q as a function of the 
inverse number of molecules N mo i for different values of the 
Coulomb interaction U. The error bars give the results of the 
quantum Monte Carlo calculations for systems with N mo i = 
4, 8, 16, 32, and, where necessary, 64 molecules. The lines 
are merely to guide the eye and identify the value of U in the 
corresponding calculations. 
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FIG. 2. Estimate of the critical ratio U c /W for a 
multi-band Hubbard model describing the doped Fm3m Ful- 
lerides (fee lattice with orientational disorder) at integer fill- 
ings n. The error bars give the estimates for U c /W from the 
quantum Monte Carlo calculations. The shaded region indi- 
cates the U /VK-range in which the real materials are believed 
to fall. 



D. Hartree-Fock calculations 
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FIG. 3. Gap E g as a function of the Hubbard interaction 
U for a cluster of 32 molecules in Hartree-Fock approximation 
for integer fillings n = 1 . . . 5. The calculations were done for 
the same cluster of 32 molecules that was used in the quan- 
tum Monte Carlo calculations. The dashed line shows the 
finite-size contribution U/N mo i to the gap (cf. eqn. (^)). The 
Mott transition occurs around U » 0.4 eV (U/W « 0.65). U c 
depends only weakly on the filling n, increasing slightly with 
increasing n. 



III. INTERPRETATION 



It is instructive to compare the results of the quantum 
Monte Carlo calculations with the predictions of Hartree- 
Fock theory. Figure |3 shows the the gap E g calculated 
for the Hamiltonian (^j) within the Hartree-Fock approx- 
imation for the different integer fillings. Compared with 
quantum Monte Carlo, the gap opens much too early, 
around U ~ 0.4 eV (U/W w 0.65). Furthermore, there 
is only a very weak doping dependence: U c somewhat 
increases with the filling — in qualitative disagreement 
with the quantum Monte Carlo results. This failure is a 
direct consequence of the mean-field approximation. In 
Hartree-Fock the only way to avoid multiple occupancies 
of the molecules, in order to reduce the Coulomb repul- 
sion, is to renormalize the on-site energy for the orbitals, 
thereby localizing the electrons in certain orbitals. For 
the Hamiltonian ([!]) this on-site energy is, apart from a 
trivial offset, given by e lma = U(J2 m '<r' n i™>cr> ~ n ima ). 
Lowering the Coulomb energy in this way will, however, 
increase the kinetic energy. For small changes in the on- 
site energies this increase will scale like the inverse of the 
density of states at the Fermi level. This suggests that 
the critical U should be the larger, the smaller the den- 
sity of states at the Fermi level. Inspecting the density 
of states N(e) for the non-interacting Hamiltonian (see 
e.g. Fig. 3 of Ref. 0), we find that this is indeed the 
case: N(e) slightly decreases with filling, explaining the 
corresponding increase in U c . Hence the weak, but qual- 
itatively wrong, doping dependence in Hartree-Fock can 
be understood as an effect of the small variation in the 
density of states of the non-interacting system. 



A. Hopping enhancement 

To find a simple interpretation for the doping depen- 
dence of the Mott transition we consider the limit of large 
Coulomb interaction U. In that limit the Coulomb en- 
ergy dominates and we can estimate the energies entering 
the gap equation (|J) by considering electron configura- 
tions in real space. According to the Hamiltonian (||) the 
contribution to the Coulomb energy from a molecule that 
is occupied by m electrons is U m(m — l)/2. Thus the 
energy of a system with filling n is minimized for con- 
figurations with exactly n electrons per molecule. The 
hopping of an electron to a neighboring molecule would 
cost the Coulomb energy U and is therefore strongly sup- 
pressed in the large- U limit. The energy for a cluster of 
N mo i molecules with N = nN mo i electrons (filling n) is 
then given by 

E(N) = n{n 2 1) N mol U + (D(t 2 /U), (6) 

where t is a typical hopping matrix element. Adding an 
extra electron increases the Coulomb energy by n U, re- 
moving an electron reduces it by (n— 1) U . But there will 
also be a kinetic contribution to the energy E(N ± 1), 
since the extra charge can hop without any additional 
cost in Coulomb energy. To estimate the kinetic en- 
ergy we calculate the matrix element for the hopping of 
the extra charge to a neighboring molecule. This ma- 
trix element will of course depend on the arrangement 
of the other N electrons. It is well known that for the 
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non-degenerate Hubbard model a ferromagnetic arrange- 
ment o£-the spins is energetically favored (Nagaoka's 
theoremtU) , allowing the extra charge to hop without dis- 
turbing the background spins. For a degenerate Hubbard 
model, however, the hopping matrix element is larger e.g. 
for an-antiferromagnetic arrangement of the background 
spinsJl3 This is illustrated in Fig. || a) for an extra elec- 
tron in a system with filling 2. Now, instead of only the 
extra electron, any one out of the three equivalent elec- 
trons can hop to the neighboring molecule. Denoting the 
state with the extra electron on molecule i by \i), we find 
that the second moment of the Hamiltonian (i\H 2 \i) is 
given by the number of hopping channels k (in the present 
case k = 3) times the number of (equivalent) nearest 
neighbors Z times the single-electron hopping matrix el- 
ement t squared. Thus by inserting the identity in the 
form where \j) denotes the state where any 

one of the electrons has hopped form molecule i to the 
neighboring molecule j, we find 



(i\H\j)=Vkt, 



(7) 



i.e. the hopping matrix element is enhanced by a factor 
of Vk over the one-particle hopping matrix element t. In 
a similar way we find for the system with an extra hole 
(Fig. [|b) a hopping enhancement of \fk with k — 2. The 
hopping enhancements for other fillings are listed in ta- 
ble [nj, where fc_ denotes the enhancement for a system 
with an extra hole, and k+ is for a system with an extra 
electron. 




FIG. 4. Illustration of how a) an extra electron or b) an 
extra hole can hop against an integer-filled background (here: 
degeneracy TV = 3, filling n = 2). For simplicitytj we consider 
the case where electrons can only hop between orbitals with 
the same quantum number, i.e. t im j m i — for m 7^ m! . 



For a single electron the kinetic energy is of the or- 
der of —W/2, where W is the one-electron band width. 
The enhancement factor \fk in the many-body case then 
suggests that the kinetic energy for the extra charge is 
correspondingly enhanced, implying 

E(N + 1) w E(N) + n U- y/k^ W/2 
E(N - 1) w E(N) -(n-l)U- yfkZ W/2. 



Combining these results we find 



E a f*U 



W. 



(8) 



i.e. the hopping enhancement leads to a reduction of the 
gap described by the factor c = (y / fc+ + y / /cT)/2. This re- 
duction is largest (ss 1.73) for n=3, and becomes smaller 
away from half-filling: c ss 1.57 for n = 2, 4, and c w 1.21 
for fillings 1 and 5. Extrapolating (|J) to intermediate U 
we find that the gap opens for U larger than U c — c W. 
Therefore the above argument predicts that the critical 
U for the Mott transition depends strongly on the filling, 
with U c being largest at half-filling and decreasing away 
from half-filling. This is qualitatively the same behavior 
as we have found in the Monte Carlo calculations. We 
note, however, that the argument we have presented is 
not exact. First, the hopping of an extra charge against 
an antiferromagnetically ordered background will leave 
behind a trace of flipped spins. Therefore the analogy 
with the one-electron case for determining the kinetic 
energy in the large- [/ limit is not exact. Second, using 
(||) for determining U c involves extrapolating the results 
obtained in the limit of large U to intermediate values of 
the Coulomb interaction, where the Mott transition takes 
place. Finally, considering only one nearest neighbor in 
the hopping argument (cf. Fig. Q) implicitly assumes that 
we are dealing with a bipartite lattice, where all nearest 
neighbors are equivalent. 



B. Origin of the asymmetry 

To analyze the asymmetry in the gaps around half- 
filling we use the following exact relation for the kinetic 
energy in the limit of infinite U , which follows from an 
electron-hole transformation 



TABLE III. Hopping enhancement for different fillings of 
a 3-fold degenerate band. 



filling 




enhancement: 


Vk 




n = 1 


fe_ 


= 1 


k+ 


= 2 


n = 2 


k- 


= 2 


k+ 


= 3 


n = 3 


k- 


= 3 


fe+ 


= 3 


n = 4 


k- 


= 3 


k + 


= 2 


n = 5 


k- 


= 2 


k+ 


= 1 



T zz (nN mol ± 1) = -T zz ({2N - n)N mol T 1). (9) 

(Note how this symmetry is reflected in the hopping en- 
hancements shown in Table III.) Since the gap for filling 
n is given by 

E g (n) = U + T min (nN mo i - 1) + T min (nN mol + 1), 

the asymmetry A = E g (n) — E g (2N — n) in the gaps can 
be written entirely in terms of energies for systems with 
an extra electron: 
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A 



-T max {{2N - n)N mol + 1) + T min (nN mo i + 1) 
-T mirl ((2N - n)N mol + 1) + T max (nN mol + 1). 



For a bipartite system the spectrum for a given filling 
will be symmetric, in particular T m i„ + T ma;E = 0, and 
thus there will be no asymmetry in the gaps: A = 0. 
Frustration breaks this symmetry. To study the effect 
of frustration we perform a Lanczos calculation in the 
large-JJ limit, starting from a configuration \v ) of the 
type shown in Fig. [|. The leading effect of frustration 
is given by the third moment, which already enters after 
the first Lanczos step. Diagonalizing the Lanczos matrix 
and expressing everything in terms of the moments .of the 
Hamiltonian, the extreme eigenvalues are given byO 



A 3 ± y/4A 



2A 2 



(10) 



where Ak = (vo\H k \vo) denotes the k th moment of H, 
and Ai = (vq\H\vo) = for a state like in Fig. ||. 
From this expression it is clear that the "band width" 
£max — £min is essentially given by the second moment, 
and that an enhancement of A 2 by a factor of k leads 
to an increase in the band width by a factor of vfe, as 
already described above. The main effect of the third mo- 
ment (i.e. of frustration) is to shift the extremal eigenval- 
ues, where the shift is determined by the third moment. 

To get a contribution to the third moment the initial 
state \vq) must be recovered after three hops. This is 
only possible if the extra charge hops around a trian- 
gle, without changing spins along its path. For a state 
with an extra electron this means that one and the same 
electron has to perform the triangular hop. Therefore, 
even in the many-body case, for each such electron we 
get the same contribution to the third moment as in the 
single electron case. It therefore makes sense to write 
the third moment A 3 (n) for a system with nN mo i + 1 
electrons in terms of the third moment A| of the sin- 
gle electron problem: A 3 (n) — fc+(n) A|, where k+(ti) 
describes the many-body effects, just like we introduced 
k+(n) to describe the many-body enhancement of the 
second moment. Using these definitions we find that the 
size of the asymmetry A in the gaps can be estimated 
by the doping dependence of the (positive) enhancement 
factors n+(n) of the third moment, while the overall sign 
is determined by the single-electron moments: 



A 



K+(n) K + (2N-n)\ A s 3 
k + (n) " k+(2N -n) ) Al, 



(11) 



To understand the doping dependence of n+/k+ we 
proceed in two steps. First we observe that the upper 
limit for the number of different electrons that can per- 
form a triangular hop is given by the number fc + of elec- 
trons that can hop to a nearest neighbor. Hence, if frus- 
tration is not suppressed, n + /k + — 1. For filling n = 1, 
N > 2 this upper limit can always be achieved without 
compromising large next-neighbor hopping by arranging 
the electrons in such a way as to avoid each other. This 



is shown in Fig. |^. For the corresponding filling 2N — 1 
the electrons can no longer be completely separated in 
that way. Thus the channel for triangular hops will be 
blocked by the the Pauli principle, reducing n + /k + . In 
that way for the larger fillings frustration is reduced. 

-f if 



n=1 

k=2,K + =2 



n=2N-1 

k=1,K=0 



FIG. 5. Typical states for a triangle with orbital degener- 
acy N = 2 and hopping only between like orbitals. For filling 
n — 1 there is no reduction of frustration (k+ = For 
filling larger than N the "background electrons" block the 
triangular moves, suppressing frustration. 

This reduction of frustration can already be seen in the 
simple model of a triangle with orbital degeneracy N = 2 
(cf. Fig. ||). Choosing matrix elements t = 1 for hopping 
between like orbitals we find for filling n = 1 a strong 
asymmetry T min (3n + 1) = -2 and T max (3n + 1) = +4, 
while at filling n' = 2N — n — 3 there is no asymmetry 
in the extremal eigenvalues: T^ax (Zn' + 1) = ±2. We 
note that flipping one spin in the configuration for filling 
n = 3 would allow for a triangular hop. In a Lanczos cal- 
culation this spin polarized configuration gives, however, 
only extremal eigenvalues T m i n = — 2 and T max = +1. 
The states described here for a triangle can be easily 
adapted to the situation in an fee lattice, where the third 
moment involves hopping to the nearest neighbor sites, 
which form connected triangles. 

From the non-interacting density of states for our 
model of the doped Fullerenes (cf. e.g. Fig. 3 of Ref. [Tl| ) 
we see that both e m i n and e rnax are shifted upwards, 
compared to the center of the band, hence, looking at 
( |l0| ) we find that for a single electron the third moment 
is positive: A 3 > 0. Together with the reduction of the 
frustration for larger filling, we therefore expect from ( |Tl| ) 
that for the alkali doped Fullerenes E g (n) > E g (2N ~n)\ 
i.e. U c (n) < U C (2N — n), as is found in the Monte Carlo 
calculations. 



IV. SUMMARY 

Using quantum Monte Carlo, we have analyzed a 
model of alkali-doped Fullerenes and found that the Mott 
transition strongly depends on the (integer) filling n. U c 
is largest for n — 3 and decreases away from half-filling. 
This result is qualitatively different from both, the re- 
sults of density functional calculations in the local density 
approximation, and the results of Hartree-Fock calcula- 
tions. The doping dependence of the Mott transition can 



G 



be understood in terms of a simple hopping argument. 
The key observation is that, due to the orbital degener- 
acy, there are more hopping channels in the many-body 
than in the single-body case, thus leading to the degen- 
eracy enhancement yk discussed above. In addition, due 
to frustration, the gaps are not symmetric around half- 
filling. 

The Gutzwiller approximation for a paranaagnetic 
state also predicts a degeneracy enhancement S3 For a 
half-filled system, the predicted enhancement is, how- 
ever, linear in the degeneracy (N + 1) instead of y/~N 
as suggested by the hopping argument, of Sec. Ill and 
as also found in infinite dimensions. E3 The results of 
the Gutzwiller approximation are reproduced by a slave- 
boson calculation in the saddle-point approximation^ In 
dynamical mean-field theory a degeneracy enhancement 
and a reduction of XL— away from half-filling, similar to 
our result, is found. [LSlIa 



ACKNOWLEDGMENTS 

This work has been supported by the Alexander-von- 
Humboldt-Stiftung under the Feodor-Lynen-Program 
and the Max-Planck-Forschungspreis, and by the Depart- 
ment of Energy, grant DEFG 02-96ER45439. 



1 N.F. Mott, Metal- Insulator Transitions, 2nd edition, Tay- 
lor and Francis, London, 1990; W.F. Brinkman and T.M. 
Rice, Phys. Rev. B 2, 4302 (1970). 

2 O. Gunnarsson, Rev. Mod. Phys. 69, 575 (1997). 

3 O. Gunnarsson, E. Koch and R.M. Martin, Phys. Rev. B 
54, R11026 (1996). 

4 T. Yildirim, L. Barbedette, J.E. Fischer, C.L. Lin, J. 
Robert, P. Petit, and T.T.M. Palstra, Phys. Rev. Lett. 77, 
167 (1996). 

5 S.C. Erwin, Buckminsterfullerenes, (Eds. W.E. Billups and 
M.A. Ciufolini), VCH Publishers, New York, pp. 217. 



6 O. Gunnarsson, S. Satpathy, O. Jepsen, and O.K. Ander- 
sen, Phys. Rev. Lett. 67, 3002 (1991); S. Satpathy, V.P. 
Antropov, O.K. Andersen, O. Jepsen, O. Gunnarsson, and 
A.I. Liechtenstein, Phys. Rev. B 46, 1773 (1992). 

7 P.W. Stephens, L. Mihaly, PL. Lee, RL. Whetten, S.M. 
Huang, R. Kaner, F. Diederich, and K. Holczer, Nature 
351, 632 (1991). 

8 I.I. Mazin, A.I. Liechtenstein, O. Gunnarsson, O.K. Ander- 
sen, V.P. Antropov, and S.E. Burkov, Phys. Rev. Lett. 26, 
4142 (1993). 

9 D.F.B. ten Haaf, H.J.M. van Bemmel, J.M.J, van Leeuwen, 
W. van Saarloos, and D.M. Ceperley, Phys. Rev. B 51, 
13039 (1995); H.J.M. van Bemmel, D.F.B. ten Haaf, W. 
van Saarloos, J.M.J, van Leeuwen, and G. An, Phys. Rev. 
Lett. 72, 2442 (1994); M. Calandra and S. Sorella, Phys. 
Rev. B 57, 11446 (1998); A.C. Cosentini, M. Capone, L. 
Guidoni, and GB. Bachelet, Phys. Rev. B 58, R14685 
(1998). 

10 E. Koch, O. Gunnarsson, and R.M. Martin, to be published 
in Phys. Rev. B (1999) BZ6420. 

11 MP. Gelfand and JP. Lu, Phys. Rev. Lett. 68, 1050 
(1992). 

12 Y. Nagaoka, Phys. Rev. 147, 392 (1966). 

13 Here we consider the case where there is only hopping be- 
tween the same orbitals on different sites. The argument 
can also be made for the case where hopping is also allowed 
between different orbitals. The resulting hopping enhance- 
ment factors \fk relative to the single-electron case are the 
same. 

14 We note that eqn. ( p^| ) seems to imply that frustration also 
slightly increases the band width. It is, however, unlikely 
that this small effect will survive in the full calculation. 

15 J. P. Lu, Phys. Rev. B, 49, 5687 (1994); J.P Lu, Int. J. 
Mod. Phys, 10, 3717 (1996); J. Biinnemann and W. We- 
ber, Phys. Rev. B 55, 4011 (1997). 

16 J.E. Han, M. Jarrell, D.L. Cox, Phys. Rev. B 58, R4199 
(1998). 

17 R. Fresard and G. Kotliar, Phys. Rev. B 56, 12909 (1997); 
H. Hasegawa, J. Phys. Soc. Jpn. 66, 1391 (1997). 

18 G. Kotliar and H. Kajueter, Phys. Rev. B 54, R14221 

(1996) ; M.J. Rozenberg, Phys. Rev. B 55, R4855 (1997); 
H. Kajueter and G. Kotliar, Int. J. Mod. Phys. B 11, 729 

(1997) . S. Kehrein, Phys. Rev. Lett. 81, 3912 (1998); R.M. 
Noak and F. Gebhard, Phys. Rev. Lett. 82, 1915 (1999). 



7 



